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In atom optics a material structure is commonly regarded as an amplitude mask for atom waves. 
However, atomic diffraction patterns formed using material gratings indicate that material structures 
also operate as phase masks. In this study a well collimated beam of sodium atoms is used to 
illuminate a silicon nitride grating with a period of 100 nm. During passage through the grating 
slots atoms acquire a phase shift due to the van der Waals interaction with the grating walls. As a 
result the relative intensities of the matter-wave diffraction peaks deviate from those expected for a 
purely absorbing grating. Thus a complex transmission function is required to explain the observed 
diffraction envelopes. An optics perspective to the theory of atomic diffraction from material gratings 
is put forth in the hopes of providing a more intuitive picture concerning the influence of the vdW 
potential. The van der Waals coefficient C3 = 2.7±0.8 meV nm 3 is determined by fitting a modified 
Fresnel optical theory to the experimental data. This value of C3 is consistent with a van der Waals 
interaction between atomic sodium and a silicon nitride surface. 

PACS numbers: 03.75.Dg, 39.20.+q 
Keywords: atom optics 



It is known that correlations of electromagnetic vac- 
uum field fluctuations over short distances can result in 
an attractive potential between atoms. For the case of 
an atom and a surface the potential takes the form 



V(r) = -% 



(1) 
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where r is the atom-surface distance and C3 is a coeffi- 
cient which describes the strength of the van der Waals 
(vdW) interaction Q . Equation^is often called the non- 
retarded vdW potential and is valid over distances shorter 
than the principle transition wavelength of the atoms in- 
volved. The significance of this interaction is becoming 
more prevalent as mechanical structures are being built 
on the nanometer scale. The vdW potential also plays 
an important part in chemistry, atomic force microscopy, 
and can be used to test quantum electrodynamic theory. 

Early experiments concerning the vdW interaction 
were based on the deflection of atomic beams from sur- 
faces. It was demonstrated that the deflection of ground 
state alkali 0] and Rydberg 0] atom beams from a gold 
surface is compatible with Eq. ^ Later measurements 
based on the Stark shift interpretation of the vdW po- 
tential were sufficiently accurate to distinguish be- 



tween the retarded V 



and non-retarded V 



forms. More recently atom optics techniques have been 
employed to measure the magnitude of the vdW coef- 
ficient C3. Various ground state [l| and excited noble 
gas atom beams have been diffracted using nano- 
fabricated transmission gratings in order to measure C3. 
The influence of the vdW potential has also been ob- 
served for large molecules in a Talbot-Lau interferometer 
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FIG. 1: A diagram of the experimental setup used. 

constructed with three gold gratings 0. 

In this article we present atomic diffraction of a ther- 
mal sodium atom beam and show that the data cannot 
be described by a purely absorbing grating. A diagram of 
the experimental apparatus is shown in Fig. ^ The su- 
personic beam of sodium atoms passes through a .5 mm 
diameter skimmer and is collimated by two 10 /jm slits 
separated by ~ 1 m. By changing the carrier gas the 
atom velocity can be adjusted from 0.6 to 3 km/s with 
^ <~ .1. The collimated atom beam is used to illuminate 
a silicon nitride grating with a period of d = 100 nm, 
thickness t = 150 ±5 nm, open width w = 50.5 ± 1.5 nm, 
and grating bar wedge angle a = 5.25±.75 degrees. All of 
the grating parameters are measured independently us- 
ing scanning electron microscope images. The diffraction 
pattern is measured by ionizing the sodium atoms with 
a hot Re wire and then counting the ions with a channel 
electron multiplier. 
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An optical description is helpful in gaining an intuitive 
picture of how the vdW interaction modifies the atomic 
diffraction pattern. To this end one should recall that 
the Schroedinger equation for a wave function -0 can be 
written as 



,9 ,, . 



V>M), 



(2) 



where m is mass, h is Planck's constant, and V is the 
potential . One can take the Fourier transform of Eq. 
[2] with respect to time and use the fact that => —iuj 
in the frequency domain to obtain 



v- 



1 



^(r,w) = 0, 



(3) 



where the dispersion relation u = - T JL has been utilized. 

L Ira 

Equationyjis usually referred to as the time independent 
Schroedinger equation. It is quite illuminating to recall 
that the Helmholtz equation [l(j for the electric field E 
is given by 



2 fc 2 )E(r, W ) =0, 



(4) 



where n is index of refraction. By inspection one can 
see that Eqs. |21and0]are formally equivalent where the 
quantities n and (l — ¥-) play analogous roles. Due 
to this fact many wave propagation methods developed 
in optics can be applied directly to matter wave propa- 
gation, being mindful of the fact that in optics ui = ck . 

While Eq. |2| can be formally solved using a Green's 
function approach, approximate solutions used in physi- 
cal optics can lead to a better understanding of how the 
vdW interaction affects atomic diffraction patterns. The 
Fresnel and Fraunhofcr approximations are commonly 
used in optics and represent a useful tool when faced with 
propagating the wave function %j) from the grating to the 
detector plane. The Fresnel or paraxial approximation is 
valid as long as the propagation distance z satisfies the 
inequality 



z > \x - £|. 



(5) 



This is certainly satisfied for our experiment since the 
diffraction angles are less than 10~ 3 radians and the or- 
ders are resolved. The Fraunhofcr or far-field approxima- 
tion goes beyond the Fresnel approximation by requiring 
that 



z > 



(6) 



where XdB is the de Brolglic wavelength of the atoms and 
tmax is the relevant extent in the aperture plane 
For the case of propagation from a uniformly illuminated 
grating of period d to the detector plane, £, m ax — > d and 
Eq. ©takes the form z > f^-. For our experimental 



setup d = 100 nm and XdB 



10 11 m, so the inequality 
However, our atom beam 



z f» 2 m ^> jgjjg m is met. 
diameter is on the order of 10 -5 m and so £ m ax — > 10~ 5 
m implying that the inequality in Eq. Hjjis not met. 

In light of the previous discussion it seems most ap- 
propriate to use the Fresnel approximation to model our 
experiment. According to the Fresnel approximation the 
wave function in the detector plane ip( x ) is related to 
that just after the grating by a scaled spatial Fourier 
transform 



(7) 
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where J^{} denotes a Fourier transform and is the 
Fourier conjugate variable to £ [Tlj . The quadratic phase 
factor in Eq. |7|accounts for the fact that the phase fronts 
have a parabolic shape before the far-field is reached. 
The wave function just after the grating is given 

by 



T(£) * comb 



(8) 



where comb (^\ is an array of delta functions with spac- 
ing d, the operator * denotes a convolution, and U (£) is 
complex function describing the atom beam amplitude in 
the plane of the grating. The transmission function of a 
single grating window T(£) in Eq. [S]is defined as 



no 



rcct 



(9) 



where rect(ar(?) = 1 when \arg\ < ^ and zero otherwise. 
The phase </>(£) accounts for the vdW interaction and its 
origin will be discussed later. This description of ipiO 
and T(£) in terms of the functions combQ and rectQ is 
standard Fourier optics notation and convenient due to 
its modular nature [llj . 

Equation |H1 can then be substituted into Eq. to 
obtain 



ip{x) oc 



oc 
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(10) 



where the summation index corresponds to the j 
diffraction order, the diffraction amplitude Aj is defined 
as 



Ai = T{n&\\ 
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and the beam profile in the detector plane is given by 



U(x) = T 



'U(0 



fr- 



(12) 
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From Eq. [TJ3 we can predict the atom intensity 

I{x) = \^{x)\\ (13) 

in the detector plane which can also be interpreted as the 
probability law for atoms. A distribution of atom veloc- 
ities can be incorporated by a weighted incoherent sum 
of the intensity pattern for each atom velocity / (x; v) 



In arriving at Eq. ^| diffraction due to abrupt changes 
in the potential V(£, z) has been ignored while the wave 
function propagates through the grating windows. This 
is a valid approximation due to the fact that XdB <C 



w 



<9£ Hoj 



mX 



P(v) oc v 3 exp 



m(v — u)' 
2k B T 



(14) 



(15) 



where the P(v) is the probability distribution function of 
velocities for a supersonic source, u is the average flow 
velocity, ks is Boltzmann's constant, and T is the longi- 
tudinal temperature of the beam in the moving frame of 
the atoms [l2|. 

One can see from Eq. QJJthat the diffraction pattern 
consists of replications of the beam shape \U(x)\ 2 shifted Taylor expanding the quantity (l 



in the region of the potential that corre- 
sponds to the diffraction orders of interest. The relation- 
ship between spatial regions of the potential V(£, z) and 
a given diffraction order will be discussed in subsequent 
paragraphs. It is also important to note that Eq. 1161 as- 
sumes that the potential V(£, z) exists only between the 
grating bars (i.e. V(£, z) = for z < — t or z > 0) and ne- 
glects the fact that the bars are not semi- infinite planes. 
Theoretical work done by Spurch et al. 01 suggests that 
the vdW potential corresponding to our ~ 50 nm grating 
bar width is very similar to that of a semi-infinite plane 
in the £ direction. Since the phase 4>(£) from Eq. 1 161 only 
depends on the integral of the potential in the z direction 
one would also expect that edge effects in V(£,z) due to 
the finite grating thickness t are a small correction. 

If the particle energy huj is much greater than the po- 
tential V(£, z) then Eq. E3 can be further simplified by 



by integer multiples of Ad g g ° with relative intensities de- 
termined by the modulus squared of Eq. 1111 An im- 
portant feature to notice in Eq. 1111 is that a diffraction 
order in the detector plane corresponds to a spatial fre- 
quency in the grating plane through the relation f^ = i- 
This highlights the connection between the spatially de- 
pendent phase </>(£) in Eq. [!|] and the magnitude of the 
diffraction orders in Eq. ITU1 

The earlier assertion that </>(£) in Eq. ^somehow incor- 
porates the vdW interaction into the optical propagation 
theory can be understood by recalling from Eq. Eflthat 

1 /2 

the index of refraction n and quantity (l — r~) play 
similar roles in optics and atom optics, respectively. In 
optics one calculates a phase shift <fi induced by a glass 
plate by multiplying the wavenumber in the material nk a 
by the thickness of the plate L (i.e. <j> = nk a L). Just as 
in the optics case one can calculate the phase shift </>(£) 
accumulated by the wave function passing through the 
grating windows 



the leading order term in 



1/2 



and keeping 



m = tk - 



V{t-,z)dz\ Tvjj-»V(£,z), (17) 



hk 

through the use of the dispersion relation lo = and 
p = mv = Hk . Equation 1171 is often called the Eikonal 
approximation. The term tk a in Eq. I17l is independent of 
£ and of no consequence in Eq. ^]so it can be neglected. 
One can see from Eq. El that if V(£,,z) — > then Eq. 
Illl rcduccs to the sine diffraction envelope expected from 
a purely absorbing grating. Furthermore, it is now clear 
from Eqs. ^2 and El that the relative heights of the 
diffraction orders are altered in a way that depends on 
V(£, z) as well as the atom beam velocity v. 

As a simple model one can represent the potenial in 
Eq. H3as the sum of the potential due to the two interior 
walls of the grating window 



( wavenumber 
yin potential 



K 1 



( differential 
Vthickness 

1/2 

dz. 



(16) 



0(0 = -T- [9-(a)V-(0 - g+(a)V + (0] , (18) 
nv 

where the function g±(a) incorporates the influence of 
the wedge angle a 



1± 



where t is the thickness of the grating and V(£, z) is the 
potential the atoms experience between the grating bars 
due to the vdW interaction. Thus the vdW interaction 
is analogous to a glass plate with a spatially dependent 
index of refraction, a kind of diverging lens that fills each 
grating window. The result in Eq. 1161 is consistent with 
the wave function phase according to the WKB approx- 
imation |9l. 



9±{ot) 



i ± 



«±¥T 



(19) 



and V±(0 = C 3 (£ ± f ) 3 is implied by Eq. □ Equa- 
tions El an d El ar e arrived at by carrying out the inte- 
gration in Eq. 1171 while assuming that the open grat- 
ing width if varies in the propagation direction z as 
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FIG. 2: Observed diffraction patterns of four different atom 
velocities. The numbers next to the peaks indicate the diffrac- 
tion order. Molecular Na2 peaks are also visible between ze- 
roth and first orders for slower velocities. 



w(z) = w + 2ztana. Since the principle transition wave- 
length of Na (590 nm) is much larger than ^ (i- e - the 
maximum atom-surface distance of ~ 25 nm) the non- 
retarded form of the vdW potential is appropriate. 

It is not immediately obvious how the phase represen- 
tation in Eq. 1 171 will affect the far-field diffraction pattern 
or if the Eikonal approximation is appropriate in light of 
Eq. HI (i.e. V±(£) -oo as |£| f ). In order to 
address this it is helpful to introduce the concept of an 
instantaneous spatial frequency 01 



I 

d' 



(20) 



where ^ is the grating window location corresponding to 
the diffraction order j as in Eq. ^2 For the limiting case 
of a — > the geometry factor g± (a) — > 1 the higher order 



terms in Eq. I16l will become important when £j — > £ c and 

C 3 (C c - §)~ 3 w fi^. If Eq. EH is inserted into Eq. EH 
with the previously mentioned limits one can solve for 
the diffraction order j c at which the approximation in 
Eq. El breaks down 



3 Aiq^ 



C 3 d 



3k a t fd 3 fu 



2 



c. 



(21) 



For the present experiment ~ 10 5 and (^^p) 3 ~ 

10s which implies that j c ~ 10 7 . Thus the approximation 
in Eq. 1171 is appropriate since we typically concerned 
with only the first ten diffraction orders. In fact, the 
paraxial approximation will become invalid before Eq. 
El becomes invalid due to the fact the diffraction order 
spacing is typically Ad ^ z ° ~ .1 mm. It is also interesting 
to note that using Eqs. 1181 andl2"Ulone can solve for the 
position in the grating window 



w ( itC^d 
2 



V jhv 



J > 1, 



(22) 



corresponding to a particular diffraction order j . If j = 1 
in Eq. 1221 then £i w 3.2 nm and since £j ~ the shape 
of the diffraction amplitude in Eq. Illl depends on a small 
region of the potential near an an atom-surface distance 
of - 20 nm. 

The experimental data for diffraction patterns of four 
different atom beam velocities are displayed in Fig. [21 
One can see from Fig. [21 that the second order diffraction 
peak is almost completely suppressed for the faster atoms 
whereas it is quite pronounced for the slower atoms. This 
velocity dependence is a clear indication that a complex 
transmission functions such as Eq. El (i.e. C3 ^ 0) is 
required to explain the data. A least-squares fit to Eqs. 
HOIandrHlis used to determine diffraction envelope \Aj\ 2 
and the average velocity. It is clear from Fig. [21 that 
the diffraction orders overlap to some extent, hence the 
tails of the beam shape are important when determining 
,4j| 2 . The broad tails of the beam shape were not ad- 
equately described by a Gaussian so an empirical shape 
using a fixed collimating geometry was derived from the 
measured raw beam profile and used for \U(x)\ 2 . 

The diffraction amplitudes \Aj\ determined from Fig. 
[21 for the various velocities are displayed in Fig. [21 The 
vdW coefficient C3 is determined by a least-squares fit 
to this reduced data with the modulus squared of Eq. 
im All of the grating parameters are determined inde- 
pendently, therefore C3 is the only free parameter. Data 
from each velocity is fit simultaneously with the same C3. 
It is clear that a purely absorbing grating (i.e. C3 = 0) 
is inconsistent with all of the observed \Aj | especially at 
lower velocities for which the phase <fi(£) is much larger. 
Uncertainty in the determination of the grating param- 
eter w and the exact shape of the potential in Eq. El 
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FIG. 3: Diffraction order intensities and best fit diffraction 
amplitude \Aj | 2 . Notice how the theory for a purely absorbing 
grating fails to describe the data. 



may be responsible for the slight deviation from theory 
evident in Fig. |3J 

A study of the systematic errors in our experiment and 
analysis suggest that w is largest source of uncertainty 
when calculating C3. One can numerically calculate the 
function Cs(w), which is the best fit C3 as a function of w, 
whose linear dependence around the physical value of w 
is found to be 5& I _ = .52 meV nm 2 . The error 

aw I w — 50.5 nm 

in C3 is arrived at by taking the product of this slope and 
the 1.5 nm uncertainty in w. After carrying out the pre- 
viously described analysis we obtain a value for the vdW 
coefficient C3 = 2.7 ± 0.8 meV nm 3 . The uncertainty 
determined this way is considerably larger than the sta- 
tistical uncertainty in C3 from the least-squares fitting 
procedure. The uncertainty due to w is also larger than 
the systematic corrections due to the atom beam pro- 
file or uncertainties due to imperfect knowledge of the 
grating parameters: d, t, and a. 



TABLE I: Measured and calculated values of C3 



Method 


Cz [meVnm 3 ] 


This experiment 


2.7 ±0.8 


Na and perfect conductor ^j] 


7.60 


Na* and perfect conductor ^3, \Q 


6.29 


Na* and Na surface 


4.1 


Na* and SiN^ surface 


3.2 


Na* and SiN^ with a 1-nm Na layer§ 


3.8 


* indicates a one-oscillator model for atomic polarizability. 
§ indicates C3 evaluated 20 nm from the first surface. 



To compare our experimental measurement with theo- 
retical predictions of the van der Waal potential strength, 
we evaluate five different theoretical cases for sodium 



atoms and various surfaces in Table [I] The Lifshitz for- 
mula 17] for C3 is 



C, 



h 

-in 



duaiiu f^l ] , (23) 

^M + i 



where a(iu) is the dynamic polarizability of the atom and 
e(iw) is the permittivity of the surface material, both of 
which arc a function of complex frequency. 

A single Lorentz oscillator model for an atom (i.e. ne- 
glecting all but the valence electron) with no damping 
gives an expression for polarizability [Isj j 



aiiuS) 



<0) 



l + (^) 2 



(24) 



For sodium atoms a(0) = 24.1 A 3 and cjq = 

2-7rc/(590 nm). Combining this with a perfect conduc- 
tor (i.e. e = 00) in Eq. 1231 gives C3 = 6.29 meV nm 3 . 
This value agrees well with the non-retarded limit cal- 
culated in reference [Hi ] for sodium atoms with a single 
valence electron. 

For more accurately modeled sodium atoms and a per- 
fect conductor, Derevianko et al. 0] calculated C3 = 
7.60 meV nm 3 and reported a range of values spanning 
0.08 meV nm 3 based on different many-body calculation 
methods which all include the effect of core electrons. It 
is noteworthy that 16% of this recommended value is due 
to the core electrons (l5| . 

For a metal surface, the Drude model describes e(iu>) 
in terms of the plasma frequency and damping: 



LU(UJ + 7) 



(25) 



For sodium metal, huj p = 5.8 eV and fry = 23 meV, 
resulting in C3 = 4.1 meV nm 3 for a sodium atom 
and a bulk sodium surface. Presumably this calcula- 
tion also under-estimates C3 because the core electrons 
are neglected. However, the calculation error is probably 
smaller than that of a perfect conductor because the core 
electron excitations are at frequencies comparable to u> p . 

For an insulating surface of silicon nitride, which is the 
diffraction grating material, Bruhl et al used a model 
with 



e(iu>) = 



lu 2 + (1 + go)u 2 
^ 2 + (1 - .9o)^o 



(26) 



where hujo = E s = 13eV and go = 0.588 is the material 
response function at zero frequency. Using Eqs. [23 EH 
and 1261 gives a value of C3 = 3.2 meV nm 3 . 

A multilayered surface makes a vdW potential that no 
longer depends exactly on r -3 , even in the non-retarded 
limit. We used Equation 4.10 from reference 0] to 
calculate V(r) for thin films of sodium on a slab of 
silicon nitride. Because our experiment is sensitive to 
atom-surface distances in the region 20 nm, we report 
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the nominal value of C3 from these calculations using 
C3 = V(20 nm) x (20 nm) 3 . Evaluated this way, isolated 
thin films make a smaller C3 as r increases. Films on a 
substrate make C3 vary from the value associated with 
the bulk film material to the value associated with the 
bulk substrate material as r increases. 

In conclusion an optics perspective to the theory of 
atomic diffraction from a material grating has been put 
forth. The results in Eqs. E]El an dEDhave been derived 
using Fourier optics techniques and appear to be consis- 
tent with the diffraction theory presented in [2(| • Diffrac- 
tion data for a sodium atom beam at four different ve- 
locities show clear evidence of atom-surface interactions 
with the silicon nitride grating. A complex transmission 
function such as that in Eq. EI is required to explain the 
data. The measured value of C3 = 2.7 ± 0.8 meV nm 3 
is limited in precision by uncertainty of the grating pa- 
rameter w. Based on the results in Table [I] for a single 
Lorentz oscillator the new measurement of C3 presented 
in this article is consistent with a vdW interaction be- 
tween atomic sodium and a silicon nitride surface. Our 
measurement is inconsistent with a perfectly conduct- 
ing surface and also a silicon nitride surface coated with 
more than one nm of bulk sodium. This implies that 
atomic diffraction from a material grating may provide 
the means to test the theory of vdW interactions with a 
multi-layered surface by using coated gratings. 

The authors would like to thank Hermann Uys for tech- 
nical assistance. 



[3 

[5 

[6 

[7; 

[§ 
[9 
[10 

[11 
[12 



[13 
[14 
[15 

[16 

[17 
[18 

[19 



[20 



[1] P. W. Milonni, The Quantum Vacuum (Academic Press, 
1994). 

[2] A. Shih and V. A. Parsegian, Phys. Rev. A 12, 835 
(1975). 



A. Anderson, S. Haroche, E. A. Hinds, J. W., and 
D. Meschede, Phys. Rev. A 37, 3594 (1988). 

C. I. Sukenik, M. G. Boshier, D. Cho, V. Sandoghdar, 
and E. A. Hinds, Phys. Rev. Lett. 70, 560 (1993). 

R. E. Grisenti, W. Schollkopf, J. P. Toennies, G. C. 
Hegerfeldt, and T. Kohler, Phys. Rev. Lett. 83, 1755 
(1999). 

R. Bruhl, P. Fouquet, R. E. Grisenti, J. P. Toennies, 
G. C. Hegerfeldt, T. Kohler, M. Stoll, and D. Walter, 
Europhys. Lett. 59, 357 (2002). 

B. Brezger, L. Hackermuller, S. Uttenthaler, 
J. Petschinka, M. Arndt, and A. Zeilinger, Phys. 
Rev. Lett. 88, 100404 (2002). 

T. A. Savas, M. L. Schattenburg, J. M. Carter, and H. I. 
Smith, J. Vac. Sci. Tech. B 14, 4167 (1996). 

D. J. Griffiths, Introduction to Quantum Mechanics 
(Prentice Half, 1995). 

J. D. Jackson, Classical Electrodynamics (John Wifey & 
Sons, 1999). 

J. W. Goodman, Introduction to Fourier Optics 
(McGraw-Hill, 1996). 

F. B. Dunning and R. G. Hulet, eds., Atomic, Molecular, 
and Optical Physics: Atoms and Molecules (Academic 
Press, 1996). 

F. Zhou and L. Spruch, Phys. Rev. A 52, 297 (1995). 
R. W. Boyd, Nonlinear Optics (Academic Press, 1992). 
A. Derevianko, W. Johnson, M. Safranova, and J. Baab, 
Phys. Ref. Lett. 82, 3589 (1999). 

M. Marinescu, A. Dalgarno, and J. Baab, Phys. Rev. A 
55, 1530 (1997). 
Lifshitz, JETP 73 (1956). 

P. Meystre and S. Sargent, Elements of quantum optics 
(1998). 

C. Ekstrom, J. Schmiedmayer, M. Chapman, T. Ham- 
mond, and D. E. Pritchard, Phys. Rev. A 51, 3883 
(1995). 

R. E. Grisenti, W. Schollkopf, J. P. Toennies, J. R. Man- 
son, T. A. Savas, and H. I. Smith, Phys. Rev. A 61, 
033608 (2000). 



